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Abstract 

Using microscopic theory, we investigate the properties of a spin current driven by magnetization 
dynamics. In the limit of smooth magnetization texture, the dominant spin current induced by 
the spin pumping effect is shown to be the diffusive spin current, i.e., the one arising from only a 
diffusion associated with spin accumulation. That is to say, there is no effective field that locally 
drives the spin current. We also investigate the conversion mechanism of the pumped spin current 
into a charge current by spin-orbit interactions, specifically the inverse spin Hall effect. We show 
that the spin-charge conversion does not always occur and that it depends strongly on the type 
of spin-orbit interaction. In a Rashba spin-orbit system, the local part of the charge current is 
proportional to the spin relaxation torque, and the local spin current, which does not arise from 
the spin accumulation, does not play any role in the conversion. In contrast, the diffusive spin 
current contributes to the diffusive charge current. Alternatively, for spin-orbit interactions arising 
from random impurities, the local charge current is proportional to the local spin current that 
constitutes only a small fraction of the total spin current. Clearly, the dominant spin current 
(diffusive spin current) is not converted into a charge current. Therefore, the nature of the spin 
current is fundamentally different depending on its origin and thus the spin transport and the 
spin-charge conversion behavior need to be discussed together along with spin current generation. 



1 



I. INTRODUCTION 



Spintronicsi^ is a novel technology that enables control of both charge and spin of elec- 
trons. To accomplish this aim, establishing methods for generation and observation of spin 
currents is an urgent issue. For generation in nonmagnetic conductors, several methods have 
been proposed. The standard way is to use the spin pumping effect in ferromagnetic-normal 
metal junctions.- - - In those systems, a spin current can be induced by the precession of the 
magnetization caused by an external alternating field and is then pumped into the normal 
metal. This spin pumping effect is widely used in experiments. A second technique is non- 
local spin injection in ferromagnetic-normal metal junctions.- 1 ^ In this case, a spin-polarized 
current is induced in a ferromagnet by applying an electric field. The spin-polarized current 
results in a spin imbalance at the interface that gives rise to a diffusive flow of spin without 
a charge current in the normal metal. A third technique is to use the spin Hall effect,— - — 
where the spin current flows in a transverse direction to an applied electric voltage in the 
presence of a spin-orbit interaction. Very recently, the spin Seebeck effect was discovered 
enabling a generation of a spin current from thermal gradients in ferro magnets. 16 Thus, spin 
current generation can be realized by several mechanisms using various magnetic, electric, 
thermoelectric, and quantum relativistic (spin-orbit) effects. 

In contrast, direct measurement of spin currents is still an open issue. Spin current 
detection has so far been performed by measuring subsidiary observables that arise from 
spin currents. The first such observation was accomplished by Kato et ai— by measuring 
optically the spin accumulation that appears as a result of spin currents at the edges of 
samples (GaAs and InGaAs) in spin Hall systems. The critical issue, however, is that 
spin currents are not conserved. Therefore, spin accumulation and spin currents are not 
in direct correspondence, in sharp contrast to charge currents that are conserved. The 
inverse spin Hall effect was proposed as a method for measuring spin currents electrically.— - — 
The idea is based on the argument that spin currents flowing in the presence of spin-orbit 
interactions induce an electric voltage through the inverse process associated with the spin 
Hall effect.— >2£>2i (The electric detection of spin transport was first demonstrated at the 
interface between ferromagnetic and paramagnetic metals by Johnson and Silsbee.-) Being 
electric, the inverse spin Hall effect has been widely utilized to detect spin currents.— 1 ^ - — 

Theoretical justification for the inverse spin Hall effect has been carried out on various sys- 
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terns.— £L2kM. Takahashi and Maekawar&^i investigated the inverse spin Hall effect because 
of nonlocal spin injection in ferromagnetic-nonmagnetic metal junctions, and demonstrated 
phenomenologically the relationship between charge current (j c ) and spin current (j a ), i.e., 
j c oc cr x j s , where cr is the spin polarization direction. Generation of the charge current 
by magnetization dynamics in magnetic junctions was studied phenomenologically by Wang 
et al. 25 and Xiao et alM By using microscopic theory, the direct connection between spin 
currents pumped by magnetization dynamics and induced charge currents was revealed in 
disordered metallic systems^>2& and in disordered Rashba systems.— ~— The result for metal- 
lic systems followed phenomenological predictions, j c oc cr x j s , whereas pumped charge cur- 
rents in Rashba systems deviated from this simple relation. The naive picture that all spin 
currents are converted into charge currents is, therefore, incorrect. What has been missing 
in this picture is the distinction between spin currents induced by some effective field and 
those arising from spin accumulation diffusion. The first one, local spin current, is a contri- 
bution proportional to a local value of the magnetization texture. The other contribution is 
a diffusive spin current which contains a long-range diffusion factor. As we will show below, 
this difference is crucial in discussing the inverse spin Hall effect but can not be addressed by 
phenomenological schemes. For applications, these two spin currents need to be considered 
separately since the direction of the local (field-driven) spin current is controlled by the field, 
while there is no way to control the direction of the diffusive spin current. 

The first aim of the present paper is to study these spin currents arising from magnetiza- 
tion dynamics (spin pumping effect) in the absence of spin-orbit interactions by providing a 
fully quantum-mechanical treatment of the conduction electrons. (We note that spin-orbit 
interactions are not essential in discussing the spin pumping effect.) Generation of spin cur- 
rents because of a precessing magnetization was first studied by Silsbee et al.- They showed 
that this precession generates an accumulation of spin at the interface between ferromag- 
netic and nonmagnetic domains and that the spin diffusion, i.e., spin current, arises from 
this spin accumulation. Tserkovnyak et alA argued that the pumped spin current can be 
expressed in the form j s = aS + bS x S, introducing phenomenological parameters a and b 
(a and b are proportional to their mixing conductances and S denotes the magnetization di- 
rection). They assumed that the spin-flip relaxation rate is sufficiently larger than their spin 
injection rate and the spin accumulation does not build up at the interface. Therefore, their 
pumping mechanism is different from that of Silsbee et al. Later, Brataas et al. considered 
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the effect of a back flow arising from the spin accumulation at the interface.- There have 
been, therefore, predicted two types of spin currents: one arising from a nonlocal diffusion 
of the spin accumulation and the other a local contribution driven by an effective field. We 
show in the present paper that the diffusive spin current originating from spin accumulation 
dominates in the case of the slow- varying magnetization (disordered system). This result is 
in agreement with the prediction by Silsbee et at— The scenario of the spin pumping effect 
without spin accumulation by Tserkovnyak et al— does not hold in the present situation, 
but it may apply to cases where rapid-varying magnetization occurs at the interface of a 
ferromagnetic- normal metal junctions. 

The second aim of the present paper is to clarify the spin-charge conversion mechanism 
based on microscopic theory and derive its conversion formula. By calculating a spin current 
in the absence of spin-orbit interactions and a charge current with spin-orbit interactions 
(treated as linear-order perturbation), we will reveal the spin-to-charge current conversion 
phenomena via the spin-orbit interaction. We will demonstrate that this conversion mecha- 
nism differs depending on the type of spin-orbit interaction, smoothness of the magnetization 
profile, and the disorder (electron mean free path). In particular, the origin of spin currents, 
namely, whether spin currents are generated by effective fields or arise from spin accumula- 
tion, turns out to be crucially important for spin-charge conversion efficiency. Depending on 
the spin-orbit interaction causing the conversion, we consider two cases: Rashba spin-orbit 
interactions and spin-orbit interaction originating from random impurity scattering in met- 
als. In the former, the dominant spin current is converted into a charge current. In contrast, 
in the latter, the correlation between the dominant spin current and pumped charge currents 
is very weak. Instead, the local spin current (a small fraction of the pumped spin current) 
is converted. Identifying the origin of spin currents is, therefore, crucial for realizing high 
conversion efficiency. 

II. SPIN PUMPING EFFECT 
A. Model 

We consider a disordered electron system coupled with localized spin S(x, t) (Fig. [Q. The 
localized spin can have spatial and temporal dependences, but we consider only the slowly 
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FIG. 1: (Color online) Schematic illustration of spin current generation. The localized spin 
S is assumed to be slowly varying in space (compared to the electron mean-free path) and time 
(compared to the electron lifetime) . Precession of the localized spin S in the ferromagnet generates 
a flow of the spin in the contiguous normal metal. 

varying case [namely, that the spatial correlation length is larger than the electron mean 
free path (see below)]. The total Hamiltonian of the system is given as H(t) = H + H ex (t), 
where 



k n k,p 

Hut) = -j^2( c i-i &c ^) • s «®- 

h.q 



(1) 



The first term H describes free electrons in the presence of spin-independent impurity 
scatterers and H ex represents the exchange interaction between the conducting electron and 
the local spin. We have introduced in momentum space annihilation (creation) operators 
cjjj for conduction electrons with kinetic energy = h 2 k 2 /2m, while u\ is the strength 
of the impurity scatterers, V is the system volume, r n represents the position of the nth 
impurity, J is the exchange coupling constant, <x characterizes a vector of Pauli matrices, 
the caret signifies a matrix, and S q denotes the Fourier transform of the localized spin (or 
magnetization). We note that impurity scattering gives rise to an elastic electron lifetime 
t. Let us stress that spin-orbit interactions are not considered in this section, since it is not 
essential to spin current generation and to spin-charge conversion. 

The electron-spin density is denned as pf(x, t) = (h/2)tr(i/j^(x,t)a a ^{x,t)) H , where 
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FIG. 2: Diagrammatic representations of the spin current density j s . Diagram (a) describes first- 
order contributions in J and (b) second-order contributions. The wavy lines denote the interaction 
with the local spin S and the gray shaded region represents the vertex correction from impurity 
scatterers. 

ip^(x) is the Fourier transform of cjP , tr denotes the trace over spin indices, and (• • 
represents the expectation value estimated for the total Hamiltonian H. The spin current 
density is defined (without spin-orbit interaction) as 

iH 3 
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where the spin current j£ has two direction components, one associated with the flow of 
spin in direction i and the other associated with a spin polarization in direction a. Here we 
used the lesser component of the path-ordered Green's function defined as G ka k , a ,(t, t') = 

{i/Wk>At'YUt)) H . 



B. Pumped spin current 

We carry out the calculation in a weak exchange coupling regime, J <C h/r. This 
assumption would be satisfied in a junction of normal metal and a ferromagnet, as shown in 
Fig. [1], since the interface is usually disordered.— The Feynman diagrams contributing to the 
spin current at first and second orders in the exchange coupling J are shown in Fig. [2j We 
assume a spatially smooth magnetization structure gl< 1, where q is the momentum of the 
magnetization and t denotes the mean free path for conduction electrons, and assume a slow 
precession of magnetization Qt <C 1, where Q is a frequency of magnetization dynamics. 
Contributions from Fig. E(a) reads 



x,t) 
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Here n; is the impurity concentration, f w = 6(ef — ftw) denotes the Fermi distribution 
function at zero temperature [9{oS) is the step function and e-p is the Fermi energy], g a (g r ) 
denotes the advanced (retarded) free Green's function given by = (g T koj )* = [hco — — 
[ih/2r)\ , and Il ra describes the diffusion ladder defining the vertex correction, 



n=0 ^ k 



The dominant first-order contribution to the spin current is calculated as 

j® a {x,t) = *^V(S°{x,t)), (5) 

where v denotes the density of states, D denotes a diffusion coefficient given as 2eFT/3m, 
and (• • •) describes long-range diffusion because of random impurity scattering 

<A(M> = U / («) 
J J q n 

where A(x, t) is an arbitrary function of both x and t. (We show details of the calculations 
in Appendix [A]) Similarly, we obtain second-order contributions in J [Fig. E](b)] as js = 
^•p(2) _|_ j| c ( 2 ) ; where the dynamic component (pumped spin current) ji is given by 

Ah 2 T 2 t 

k,k',q,q' w,f2,Q' 



X 



ImekglJgU 2 
2uJ 2 r 2 D 



r Vl^ilx^i)]"), (7) 
and the equilibrium component (spin super- current) jT^ is given by 

k,q,q' uj,Q,Q' 

This equilibrium flow is a spin super-current and arises from the angular difference between 
two localized spins (or magnetizations) S\ x S2, i.e., from the spin Josephson effect.— In 
Eqs. ([5]) and ([7]), the dynamic component without vertex correction does not contribute 
because it is smaller than that with vertex correction by 1 order of magnitude (q£) 2 <C 1. 
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As a result, the pumped spin current for a smooth- varying magnetization texture is a long- 
range diffusive flow (Eqs. (jSJ) and ©). From Eqs. (jSJ) and ([7]), the pumped spin current, 
js = ji 1 ^ + Js , can be written as a gradient of an effective potential ji 



Si 



j(p) a (x,t) = -V/i«( : r,t), (9) 

where 

OM) - -^^<^(a,,t)> + ^^<[S0M) x $(*,*)]«) (fO) 

This effective spin potential arises from dissipations of magnetization energy S x S and S. 
Our result, Eqs. (jUJ) and (ITU]), indicates that the sources of the spin current are S and S x S 
and this result appears to agree with phenomenological predictions for the spin pumping 
effect. 4 However, the spin current here has been diffusion-averaged and is not simply a local 
function of the source. 

Similarly, spin density is given by 

tfOM) = -^^V 2 (^0M)> + *^([ S (x,t) x S(x,t)] a ). (11) 

It satisfies the spin diffusion equation: 

p:( X ,t)-Dv 2 P :(x,t) = r«(x,t), (i2) 

where t p represents the spin relaxation given as 

T°(x,t) ee -^V 2 S a (x,t) + ^L[ S (x t t) x S(x,t)] a . (13) 

The effective spin potential is written as 

/!?(*, t) = Dtf - ^S-fof). (14) 

Specifically, it is the spin density excluding the direct spin polarization because of the ex- 
change interaction, p s — (hvJ/V)S, multiplied by D. Therefore, the total spin current is 
reducible to 

h 2 vJ 2 



12me F V 



[S(x, t) x VS{x, t)] a - Vfif(x, t). (15) 
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The last term, a diffusive contribution, can be written in terms of the spin density by using 
Eq. (j!4p and we finally obtain 

j?(x,t)=jW a (x,t)-DV P Z(x,t), (16) 

where we have defined the total equilibrium spin current as 

hu IT) h 2 u T 2 

j^"(x,t) = ^VS"(x,t) x VS{x,t)] a . (17) 

As seen in Eq. ( jT6]) . the dominant contribution of the dynamic spin current derives com- 
pletely from a diffusion of the spin polarization. Equations (I16p and (TTTj) are the main results 
describing the spin pumping phenomena. 

For a charge current j c , it is described generally as 

j c (x, t) = a c E(x, t) + j BC (x, t) - BVp c (x, t), (18) 

where a c is an electrical conductivity, E represents an electric field, j sc is a super-current 
which exists in an equilibrium situation (corresponding to jjf ^ of the spin current), and p c 
denotes a charge density. Comparing the two expressions, Eqs. ([IB"]) and ( 118"]) . we imme- 
diately see a striking difference between the charge and the spin, namely, there is no field 
that induces a spin current at least in the present perturbative regime without spin-orbit 
interactions. In other words, the spin pumping effect does not directly generate the spin 
current itself, but causes a spin imbalance or spin accumulation that then gives rise to a 
diffusive spin current. This result is in agreement with the study by Silsbee et al.,— while 
the prediction by Tserkovnyak et al.,— stating that spin pumping is not associated with spin 
accumulation, does not hold in the smooth spin case. The diffusive spin current represented 
by the vertex correction found here would correspond to the phenomenologically discussed 
back-flow because of spin accumulation at the interface.- 

Absence of an effective field for the spin current is crucial in spintronics. In fact, in charge 
electronics, the charge and its current are independently controllable and can be measured 
by different mechanisms, for instance, by capacitance means for charges and Ampere's law 
for currents. This is not the case in spin transport phenomena. If we use the spin pumping 
mechanism in a disordered system, the spin current is always accompanied by spin accumu- 
lation according to Eq. ( IT5"|) . 
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Fortunately, we know that there is an effective field for a spin current if we use spin-orbit 
interactions, specifically, as in the spin Hall system. By including spin-orbit interactions, 
the spin current is generalizable to 

%(x,t) = (TsaeifrEjfaQ+j^fat) - BV lP f(x,t), (19) 

where <tsh represents the spin Hall conductivity proportional to the spin-orbit interaction 
and eij a denotes the Levi-Civita antisymmetric tensor. Therefore, charge and spin currents, 
Eqs. ( I18p and (I19p . now look symmetric. However, there appears a crucial difference when 
one includes spin-orbit interactions, namely, the violation of the conservation law for spin. 
In fact, spin and its current in the presence of spin-orbit interactions satisfy the identity 

p:( X ,t) + v-j:(x,t) = % a (x,t), (20) 

where % represents the spin-relaxation torque arising from the spin-orbit interaction.— 
(This equation is equivalent to the diffusion equation for the spin density.— - —) The non- 
conservation of spin causes definitional ambiguity of a spin current. Since a definition of spin 
current is absolutely related to a definition of spin relaxation torque, spin currents cannot be 
defined uniquely under the spin-orbit interaction. (As a possible solution, a gauge covariant 
derivative was proposed. 3 -) Because of the spin relaxation torque %, the measurement of 
the spin current can never be carried out by simply measuring the spin density induced at 
the edge since the induced spin current can disappear because of the term % while being 



cannot be con- 



transported. (In this sense, the observation of the spin Hall effect in Ref. 
sidered as a direct observation of the spin current.) At present, there has been no indication 
of an emission of an observable field (either electric or magnetic) from the spin current and, 
therefore, in contrast to charge currents which are observable by detecting an Ampere's field, 
electromagnetic detection of spin currents is not possible. [The absence of electromagnetic 
fields induced by spin currents is reasonable from Maxwell's equations because the equations 
as determined by special relativity and U(l) gauge invariance cannot be modified by the 
spin current.] Here a serious dilemma for spintronics arises. Specifically, spin currents and 
spin densities are independently controllable only if one switches on spin-orbit interactions, 
but such interactions make it impossible to detect the spin current by measuring the spin 
accumulation. 
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FIG. 3: (Color online) Relationship between pumped spin current and interface condition in 
a ferromagnetic-normal metal junction, (a) The diffusive spin current is dominant where the 
magnetization structure S varies slowly as compared to the electron mean free path £ in normal 
metal, (b) The local spin current is dominant where the magnetization structure varies rapidly at 
the interface. 

C. Diffusive spin current vs. local spin current 

In the spin pumping effect, we have demonstrated that the diffusive spin current is dom- 
inant in a slow-varying magnetization structure, subject to £ <C A (£ is the electron mean 
free path and A represents the length scale of the magnetization structure) as depicted in 
Fig. |3](a). In an actual spin pumping system, this condition would be satisfied since the 
mean-free path in such experimental systems is known to be very short, for example, less 
than 1 nm size in Ni 81 Feig/Pt sample.— Strictly speaking, a local spin current can also exist 
but this contribution is negligibly small compared with the diffusive spin current. There- 
fore, in this situation, it is impossible that local and nonlocal spin currents coexist, as found 
in a study by Brataas et al.— The situation can be different with a clean sample or very 
sharp magnetization profile, subject to £ ^> A, as shown in Fig. |3](b). In this case, the 
local spin current without spin accumulation could become dominant and the diffusive spin 
current arising from spin accumulation is small. We therefore expect that the mechanism of 
Tserkovnayk et al. without spin accumulation may be valid. The spin current induced by 
the spin pumping effect, therefore, depends much on the disorder or the electron mean-free 
path. 

In this section, we did not take into account spin-orbit interactions. It is essential when 
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we need to discuss spin-charge conversions, which we do in the following section. 



III. SPIN-CHARGE CONVERSION 

We now consider the conversion mechanism of the pumped spin current into the charge 
current through the inverse spin Hall effect. In this section, we consider two types of spin- 
orbit interactions, the Rashba type and such interactions because of random impurity scat- 
tering, and we finally obtain the exact spin-charge conversion formula. 

A. Rashba spin-orbit interaction systems 

We first consider charge currents driven by Rashba spin-orbit interactions. The Rashba 
spin-orbit coupling was first found as a peculiar effect in a two-dimensional electron-gas sys- 
tem.— However recent studies have revealed that the Rashba system appears quite generally 
at surfaces of nonmagnetic materials without inversion symmetry and that this Rashba 
coupling can be quite large.—"— Therefore, there is a possibility that these surface effects 
contribute greatly to the inverse spin Hall effect in ferromagnetic-normal metal junction 
systems. 

The microscopic theory in the presence of Rashba spin-orbit interactions was demon- 
strated by Ohe et al.— They considered a two-dimensional electron-gas system, where the 
maximum number of diffusion ladder needs to be taken into account. The result was ex- 
tended to a three-dimensional and spatially dependent Rashba coupling case.— However, 
the account payed little attention to the charge conservation law because the analyses were 
intended only to see whether the spin-charge current conversion indeed occurs or not. In the 
following, we consider a three-dimensional system and evaluate the dominant contribution 
including one diffusion ladder as a vertex correction. By considering the vertex correction, 
we maintain charge conservation throughout the calculation. 

The Rashba spin-orbit interaction is given by 



k 

where E so describes the spin-orbit field (or strength of the Rashba coupling). In the presence 
of this interaction, the anomalous velocity resulting from the spin-orbit interaction modifies 
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FIG. 4: Dominant contributions to the pumped charge current by magnetization dynamics in a 
Rashba system: (a) is the normal charge current, j°, and (b) describes corrections to the charge 
current arising from the Rashba coupling, Double lines represent Rashba spin-orbit interactions 
(SOIs). The Rashba spin-orbit coupling gives the anomalous velocity to the conduction electrons, 
thereby modifying the definition of the charge current. 

the charge current. By denning the charge density as p c (x, t) = —etr(il)>(x,t)i/)(x,t)} H , the 
charge current density is given by j c = j° where the normal charge current is 

ieh 2 
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and the correction from the Rashba coupling is defined as 
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(23) 



We treat the Rashba spin-orbit interaction perturbatively by imposing the constraint 
E so kp ^ h/r, with k-p being the Fermi wavelength. For slow- varying magnetization textures, 
the contribution from both the first-order Rashba and the first-order exchange interactions 
vanishes identically.— 1 ^ (Strictly speaking, a contribution proportional to VE so VS arises if 
Rashba spin-orbit interactions are inhomogeneous.— ) We now consider the charge current 
to first order in the Rashba coupling and second order in the exchange coupling. 

In Fig. HI we present the Feynman diagrams for the dominant contribution which is 
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calculated (see details in Appendix [B]) as 

k,k',q,q' u),Q,n' J 



4ez/ J 2 t 2 r • ~i 

{#so x [5(a;,t) x 

ffly V '{ V '< Eso>< xS(x,t)]>). (24) 



This pumped charge current is expressible in terms of the dynamic component of the spin 
relaxation torque 7^ and the charge density^ 

Jci (x, t) = e ija afT^ )a (x, t) - DV iPc (x, t), (25) 

where a R = —2erE so /h 2 and 



T B ^(x,t) = ^[S(x,t)xS(x,t)] a , 
2uJ 2 t 2 

p c (x, t) = — —a R ■ < V x [S(x, t) x S{x, t)]). 



(26) 



We first note that Eqs. (1251) and (126]) indicate that only the damping of the local spin, 
S x S, is converted into a charge current by Rashba coupling. The equilibrium spin current, 
Eq. (|17p . does not contribute, as is reasonable from energy conservation considerations. 
Since j c is expressible in terms of %^ dyS> and p c , the naive formula for the inverse spin Hall 
effect, i.e. j c proportional locally to j s , does not hold in disordered Rashba systems. Instead, 
Eq. ( 1251) indicates that the generation mechanism of the local charge current in a Rashba 
system is the inverse effect of the spin-relaxation torque. 

The non-local part in Eq. (1231) arises from a diffusion of the charge density, which is 
written in terms of the dynamic component of the pumped spin current [Eq. (17j)] ; 

p c (x,t) = ~e ija a?jf )a (x,t). (27) 

Therefore, the (diffusive) spin current induces a charge polarization via Rashba spin-orbit 
interactions, but not a charge current itself. Since diffusive spin currents dominate in slow- 
varying magnetization structures, we expect that the ratio of the charge current to the 
pumped spin current is high [see Fig. 0a)]. 
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FIG. 5: (Color online) Conversion mechanism of the pumped spin current into a charge current via 
spin-orbit interactions (SOIs). (a) In the Rashba system, the diffusive spin current is just converted 
into a charge current, (b) With spin-orbit interactions caused by random impurity scattering, the 
diffusive spin current is present but not converted. The local spin current is converted into a charge 
current. 

B. Random impurity-induced spin-orbit interaction systems 

We now focus on spin-orbit interactions caused by random impurities. This interaction 
is denned as 



iui\ s 
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where A so is the spin-orbit strength. We have here assumed that the random impurity spin- 
orbit interaction arises from the same impurities giving rise to the electron lifetime r. In 
this case, the correction current resulting from the spin-orbit interaction is given by 

eui\ s 



f c °(x,t) 

As shown in Ref. 
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, the charge current is derived in the form 
j ci {x,t) = a im % Q £ ocal)a (x,i) - DV iPc (x,t), 



(29) 



(30) 



where a imp = 2e\ so kp/eFT and 
HuJtD 
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e ija ([ViS(x,t) x VjS(x,t)] a ). 



(31) 
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We have here denoted the local spin current as ji , . Equation (1301) seems consistent with 
the naive formula of the inverse spin Hall effect, j c oc a x j s . However, one should note 
that this local spin current is a very small correction to the dominant spin current [given by 
— DVp s in Eq. (1161) ]. Therefore, most of the spin currents generated by the spin pumping 
effect are not converted into a charge current; only a small fraction [of order (q£) 2 <C 1] 
develops into a charge current [Fig 0(b)]. Thus, the spin-charge conversion efficiency is 
small in the presence of random impurity-induced spin-orbit interactions. 

C. Conversion mechanism 

By comparing Eqs. fl25|) and (130]) to the general formulation of the charge current 
[Eq. (fT8|) ]. the spin relaxation torque 7^ dy ^ and the local spin current js° ca ^ act as effective 
electric fields for Rashba and random impurity-induced spin-orbit interactions, respectively. 
Therefore, as we mentioned above, the spin current changes only the constitutive relations 
associated with electromagnetic fields but does not change the Maxwell's equations them- 
selves. Correctly, there is a spin current caused by spin-orbit interactions. This spin current, 
however, produces a second-order charge current with respect to spin-orbit interactions at 
least. It should be negligible compared to the above results. 

From the above results, we see that spin accumulation at the interface plays a crucial 
role in determining spin pumping and spin-charge conversion mechanisms. In fact, the 
pumped charge current is proportional to the spin accumulation when Rashba interactions 
are present, but does not occur with the spin accumulation for random impurity- induced 
spin-orbit interactions. Therefore, measuring spin accumulation at the interface would pro- 
vide impetus to determine the dominant spin-orbit interaction and to clarify the spin-charge 
conversion mechanism. 

IV. CONCLUSIONS 

We have studied aspects of spin pumping and the spin-charge conversion mechanism 
through spin-orbit interactions in the disordered electron system. We showed that the spin 
current generated by the spin dynamics is a diffusive process arising from a dynamics- 
induced spin accumulation. There is, therefore, no effective field that drives the spin current 
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directly in the disordered case. We have confirmed that a charge current is induced by 
these spin-orbit interactions. This process involves the conversion of a pumped spin current 
into a charge transport, but the mechanism has turned out not to be as simple as an 
earlier phenomenological proposal, j c oc (cr x j s )iL2JL2i had anticipated. In fact, the spin- 
charge conversion depends largely on the type of spin-orbit interaction. For Rashba spin- 
orbit interactions, the charge current is given by a local contribution proportional to the 
spin relaxation torque and a diffusive contribution arising from the diffusive spin current. 
Therefore, the naive formula for the inverse spin Hall effect does not hold in the Rashba 
systems. In contrast, for random impurity- induced spin-orbit interactions, the local part 
of the charge current is written as a very small fraction of the spin current [smaller by 
0(£ 2 /X 2 ), where I and A are the electron mean-free path and the coherence length scale of 
the magnetization, respectively]. However, the dominant spin current is not converted into 
a charge current in the presence of impurities. Thus, the naive inverse spin Hall effect does 
not occur either. Our result indicates that the spin-charge conversion formula as proposed 
earlier using phenomenological arguments is too simple and the whole phenomenon needs 
discussing together with the origin of spin currents. 



Appendix A: details of calculations for the spin pumping effect 



We perform here calculations of the pumped spin current using standard perturbation 
expansion techniques. We treat the exchange coupling up to the second order. The electron 
spin and its relaxation torque densities are defined in terms of Green's function as 

Hi 



2V 
ihJ 



tr 



(M) 



k,q 



T s a (x,t) = -^£ e -^tr{[<T x S(x,t)] a d^_ iMi (t,t)}, 

k,q 



(Al) 



respectively. Before carrying out the calculation, we introduce the Dyson equation: 

Gka,k'a' (t, t) =S kk '5 aa 'g ka (t, t) 

+ V I dT J2T, eiprn ^ a (t,T)G k+paikla/ (T,t') 

^ C K n p 

-J dTj2J2^a{t,T)[a ar S q {T)]G k+q ^ a/ {T,t'), (A2) 



P <7 
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r n a 




FIG. 6: The lesser component of the Green's function involving the vertex correction can be 
divided into the three terms. The diagram is partitioned off by defining the diffusion ladder as 
the boundary. We denote contributions from the left-hand side as T, the middle representing the 
diffusion ladder as II, and the right-hand side as A. The spiral line represents all interactions and 
therefore the contribution from A depends on the interactions. 

where G kct)k > a >(t,t') = -(i/H)(T CK [cfc a (*)4'a'(*')] ) H ( T c K bein g the path-ordering operator 
denned on the Keldysh contour Ck) and g is the free Green's function which is obtained 
from the free Hamiltonian H Q . The Dyson equation is very useful in carrying out the per- 
turbation expansion because this equation can be solved iteratively. Here we assume a weak 
exchange coupling regime, J <C ft/r, and therefore we can treat the exchange interaction 
perturbatively. To evaluate the lesser component of G(t,t') = f c dTGi(t,T)G2(T,t'), we 
use the following 44 




dT Gf T) (t, T)Gf ] (T, t') + G? T) (t, T)Gf ] (T, t') . 



(t,T)Gl(T,t') , 



(A3) 



The lesser component of the free Green's function satisfies g£ u = f w {g%u — sfLj)- 
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1. First-order calculations in J 



First, we show the calculation of the spin current to first-order in the exchange interaction. 
The diagram in Fig. |5^a) is written as 



ji 1)a (x,t) 



ih 3 J 



2mV 



E 

n=0 



V 



2\ n 



<? {M?= 



I [ / . dt i+1 g ki _^(ti,t 



i=0 " b K 



i+l, 



[<T • S q (t 



n+1 



n 



i=0 " °K 



dkn+i+l9k n -i+% (tn+i+l, tn+i+2) 



(A4) 



*0=*2n+2=* 

By using tr(<j Q (j /3 ) = 25 a ^ and taking the lesser component, the equation reads 
iff J 



EV^ c -iq-x+iQt na . / a aa(l) A rr (l) 1 A ra (l" 

k,k' : q uj,Q 

I paa rraa A aa (l) prr TT rr A rr ^ _l_ P ra TT ra A ral 

"r 1 fe,qr;w,n ii <j;aj,f2 Yi -fc',q;w,n 1 fc,<2;w,C ii qr;a;,f2 iV fc',q;w,n 1 fc,q;w,C ii <5f;w,n Yi -fc'. 



(AS) 



A diagram involving the vertex correction is divided into three parts: the left-hand side, 
the diffusion ladder (middle) , and the right-hand side of the diagram shown in Fig. [6j The 
contribution from the left hand side is given as 



racial 



(rr) Tliltj a ( r ) a ( r ) 



A'. </:_'. < ..' 



pra _ n i u j _ r 



(A6) 



The diffusion ladder arising from the vertex correction is written as 



n aa(rr) 



rrra 
L1 q;u),n 



Wy; r aa(rr) 

n=0 ^ fe 

00 ✓ 



fc,<j;o;,n i ' 



(A7) 



n=0 v fc 

The contribution from the right-hand side of a diagram depends on the diagram and in 
Fig. [2(a) is given by 



. aa(rr)(l) r a(r) a(r) 



2 ' 2 1 2 ' w 1 2 



(A8) 
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Assuming slow dynamics Qt <C 1 and a spatially smooth local spin structure < 1, we 
obtain the result 



2K A J 
3mV 

huJrD 
V 



k,k',q lu,Q 



V 



Here we note simplifications in the k summation 



SL\2 



2tivt 

IT' 



(A9) 

(A10) 
(All) 



where we have put gk = gk,frw=e F - The diffusion ladder arising from the vertex correction is 
also given as 



n 



aa(rr) 



LL q,0,Q - 1 T/ 



(A12) 



n=0 



fc 



(l-^)^-^Im^(^) 2 



3m 



n=0 



(A13) 



Dq 2 r + ittr' 

Since the product of only g a (or g r ) gives a very small contribution which is of order 1/ef 
compared to the coefficient of g a and g r , the diffusion ladder reduces approximately to 
unity. We cannot, however, ignore the product of only g & (or g r ) completely because that 
contribution corresponds to an equilibrium current. 

Similarly, the spin density is also calculated in the form 
iH 2 J 



pi 1)a (x,t) 



V 

ih 2 J 



\^ \^ -iq x+iQt qa / »aa(l) .rr(l) _i_ TT ra A 1 ^ 1 ) ^ 

/ j / j e °q,n \ Iv k,q;u),Q Iy k,q;iv,Q ' LX q;u3,Sl Iv k,q;u),Q J 



k,q U),Q 



e%q x+%Ut Sq,nfuj ( I" n j gl„ui9k, 



V 

k.q U3,\l 



V 



V 2 (S a (x,t)). 



(A14) 



To first order in the exchange coupling, the spin-relaxation torque corresponding to diagram 
Fig. [7(a) vanishes as a consequence of tr<r a = 0. Therefore, the pumped spin current to first 
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(a) (b) 




FIG. 7: Diagrammatic representations of the spin-relaxation torque of the conduction electrons, 
(a) This contribution comes from first-order terms in J but vanishes, (b) The leading contribution 
arises from second-order exchange interaction terms. This contribution corresponds to local spin 
damping S x S. 

order in J follows 

p^ a (x,t) + V-j^ a (x,t) = 0. (A15) 
Hence, the spin of conduction electrons is conserved. 



2. Second-order calculations in J 

Here we derive the spin current to second order in the exchange coupling as shown in 
Fig. [2(b). This is calculated in the same manner as the first-order calculations. The pumped 
spin current reads 



fc3 t2 

k,k',q,q' CJ,n,Cl' 

A rr (2) .ra(2) i F ra TT ra A ra ( 2 ) 

ly k,q,q';u,Cl,W ' ^k,q,q';u,Q,Q' ' 1 k,q;uj,n 11 q;uj,n ly k',q,q';u},n,a' 



h 2 J 2 



3mV 



h,h',q,q' to,Cl,Cl' ^ 



4rfl' \lme k g T k Jgl 



a \2 



/r// /2 -[S(x,t) x VS(x,t)] a - 2uj2 J 2B V([S(x,t) x S(x,t)] a ). (A16) 



12me F V v /J 1/ 

In this case, we should pay particular attention to tr(a a a^a" f ) = i2e Q( g 7 and the contribution 

from the right-hand side of the diagram being given as 

i aa(rr)(2) p a(r) a(r) a(r) 

A k,q,q>;u,U,n> =Ju-(+)%9 k _i,.,_q9 k g-V^. n-an' fffc + a„, + S' 
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The pumped spin current contains an equilibrium flow given by estimating (g a ) 2 , 



£^) 2 - (A18) 



The spin density is calculated in a similar manner 
V 



k,q,q' lu,Q,Q' 



/.aa(2) __ Arr(2) _|_ TT ra A ra (^) 

A I - (1 -Jfc 1 g 1 g';w ) n,n' IY k,q,q';uj,n,n> ' LL q;bJ,n lv k,q,q';uj,Q.,W 



k,q,q' 



IvJ 



2 r 2 



r x [S(x,t)xS(x,t)] a ). (A19) 

To second-order in the exchange coupling, the spin of conduction electrons is not conserved 
and follows the general spin continuity equation 



(2)«, 



x, t) + V • j s (2)Q (a;, t) = 7; (2)a (a;, t). (A20) 



Here the spin relaxation torque is represented by the diagram in Fig. [3(b) and we obtain 
i2HJ 2 
~V 



fc,qr a;,Q 

6m 



fc,q a),Q 



- ~ rSS 7 [5(a; ' t} x V25(;r ' t)] a + [S{x ' t] x t)] °" (A21) 

Appendix B: Magnetization-pumped charge current in a Rashba system 

Here, we calculate the charge current stemming from magnetization pumping in the 
presence of Rashba spin-orbit interactions. In the Rashba system, the Dyson equation is 
modified as 

(t,t') =S kk >5 

aa' 9ka{t, t') 

+ V I dT J2J2 etP ' r "9ka(t,T)G k+pa , kW (T,f) 

n p 

— J dT'y^ j 'y^ j g ka {t,T) \(T a p ■ S q (T)\Gk+qP,k'a'(T,t'). 

^ C K p q 

- dTj2J29ka(t,T)[E so -(kx(T a ^]G k ^ kW (T,t'). (Bl) 

q 
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Since the spin-orbit interactions are generally weak, subject to E so kp <C h/r, we perform 
the perturbation expansion with respect to both the exchange interaction and the Rashba 
spin-orbit interaction. By iteration, we obtain the leading contribution shown in Fig. HJ 

2eJ 2 



jci (•") tj 



v 

h 2 ki 



E E 

k,k',q,q' u},Q,n' 
i ( 7 aa(2) 



rn 



-iq-x+iQ,t 



Aaavzj Arr(2) 7 ra(2) 

ly j]k,q,q';w,Q.,W iv j;ft,g,g';u>,n,f}' ' ly j;k,q,q';u,n,Cl 



E so X {Sg'pi X S'q-q'^-n/) 

J J 

i pra TTra A ra ( 2 ) 

' 1 k,q;uj,n 11 q;uj,n- lY j;k',q,q';u},n,n' 



if / A aa(2) .rr(2) .ra(2) 

U *i ^^fe.q.q'itj.n.n' i4 -fc,g,g';o;,Q,n' ~T IV k,q,q' 



ieJ 2 r 

i 

mV 



E E' 

k,k',q,q' w,Cl,W 



-iq-x+iQt 



E so x (5q' 5 n' x S q -q' t n-n>) £1' f^g ktUJ gl 



x 



i-rra 



|^; so x [5(aj,t) x flf(aj,t 



4eu 7 2 r 3 D r 

^ V,{V-(g so x [8(x,t) x S(x,t)]) 

(B2) 



Here we have defined contributions from the right-hand side of the diagrams as 



A 



aa(rr)(2) 

i;k,q,q';LO,n,Q' = Ju)-(+) 



mj)_/ a ( r ) a ( r ) 



a(r) 
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a(r) \ a(r) 



7ra(2) 
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= if n-2n> - 



m d 



(B3) 



+ 



m (9 

^2 9fc ,w-§5 , fe-2^ ja ,_£^^fe+2, w+ § 



2 ' 2 
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The charge density is written in terms of the lesser component of the nonequilibrium Green's 
function 



PcOM) = yE e_i? ' Itr ^ 2 , fc+f (M). (B4) 



k,q 
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In a similar manner to the charge current, the charge density is calculated as 
2eHJ 2 



Pc(x,t) 



V 



k,q,q' ui,Q,fl' 



( laa(2) _ Trr(2) 7ra(2) 

A I 1 H;k,q,q';ui,U,Q' IY i;k,q,q';u!,n,n' ' 11 q;uj,n 1 ^i;k,q,q';ui,n,n' 

4eJ 2 r 2 



E E - 

k,q,q' u),Q,£l' 



-iq-x+iQt 



E so x (-*v if y x Sq-q^n-w) q&' fl~R r q- w ^9k,w9k,<, 



4evJ 2 T 3 
h 2 V 



V-(E so x [5(aj,t) x S(x,t)]). 



(B5) 



The charge and its current densities that we have obtained satisfy the following charge 
conservation law: 

p c (x,t) + V-j c (x,t) = 0, (B6) 
indicating that our calculation has been performed correctly. 
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